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This paper is devoted to the study of two-point correlation function of the energy-momentum 
tensor {T12T12} for 5C/(2)-gluodynamics within lattice simulation of QCD. Using multilevel algo¬ 
rithm we carried out the measurement of the correlation function at the temperature T/Tc — 1.2. 

It is shown that lattice data can be described by spectral functions which interpolate between hy¬ 
drodynamics at low frequencies and asymptotic freedom at high frequencies. The results of the 
study of spectral functions allowed us to estimate the ratio of shear viscosity to the entropy density 
ri/s = 0.134 ±0.057. 
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Introduction 

One of the most important result obtained at RHIC experiment is the measurement of the elliptic flow of final 
particles [U [2]. The value of the elliptic flow measured at RHIC can be explained within the approach based on 
hydrodynamics [31], if one assumes that quark-gluon plasma (QGP) obtained after the collision is almost superfluid. 
In particular, numerical simulations of the relativistic liquid showed [^ that the upper limit on the ratio of shear 
viscosity to entropy density is rj/s < 0.4. Preferred range for the ratio is rj/s = (1 O 3) x l/ 47 r, which is very close 
to the result of N = 4 Super Yang Mills (SYM) theory at strong coupling r]/s = l/ 47 r| 7 ]. From these facts one can 
conclude that theoretical prediction of the ratio ry/s is a very interesting and important problem. 

It is known that QGP is strongly interacting system and today there are no analytical approaches which allow to 
study such systems without additional assumptions. For this reason one of the main approaches which can be used 
to study the properties of QGP and which is based on the first principles is lattice simulation of QGD. 

Lattice simulation of QGD aimed at the calculation of shear viscosity of the S'C/(3)-gluodynamics was carried out 
in papers [SHIT]. Despite rather large uncertainties one can state that the ratios 77 /s obtained in papers [101 [H] are 
close to the N = 4 SYM prediction 77 /s ~ l/ 47 r. The calculation of shear viscosity in QGD with dynamical quarks is 
still a challenging problem. 

The closeness of the SU(3)-gluodynamics ratio 77 /s to the N = 4 SYM prediction allows us to ask the question 
whether this property is the property of the SU(3)-gluodynamics or it is common to all non-abelian gauge theories. 
To address this question in this paper we are going to perform the calculation of shear viscosity of SU (2)-gluodynamics 
within lattice simulation. 

We have already performed the first measurement of shear viscosity of the S't7(2)-gluodynamics in jT^j. In this 
paper we extend our investigation to larger lattice size and apply different method to determine numerical value of 
the ratio 77 /s. 
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Details of the calculation 


Shear viscosity is related to the Euclidean correlation function of the energy-momentum tensor = \S)iuF^pF^p — 
(here we omitted for simplicity trace anomaly): 


C{xo) = T-^ J d3x(ri2(0)Ti2(xo,x)), 


( 1 ) 


where T is the temperature of the system. The correlation function 0 can be written in terms of the spectral function 
p{(jj) as follows 


C{xo) = T 
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The spectral function contains a lot of important information about the properties of medium. In particular, to 
find shear viscosity from spectral function one uses the Kubo formula 
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Lattice calculation of shear viscosity can be divided into two parts. The first part is the measurement of the correlation 
function C{xo) with sufficient accuracy. This part of the calculation requires large computational resources but for 
the gluodynamics the accuracy of the correlator can be dramatically improved with the help of the two-level algorithm 
M- The second part is the determination of the spectral function p{uj) from the correlation function C{xo). The 
last part of the calculation is probably the most complicated, since one should determine continuous spectral function 
p{u}) from integral equation ([^ for the set of 0(10) values of the function C{xo) measured in lattice simulation. 

To approach the solution of the integral equation 0 one should take into account the properties of the spectral 
function. Important properties of the spectral function are positivity: p(ijj)/uj > 0 and oddness: p{—io) = —p{uj). It is 
also important to write the expression for the spectral function at the leading-order approximation in strong coupling 
constant [15] 
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where cIa = — 1 = 3 for the SU(2)-gluodynamics. One also knows next-to-leading order expression for the spectral 

function at large uj |16j 
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It should be noted here that at large oj the spectral function scales as p{uj) ^ what leads to a large perturbative 
contribution to the correlation function for all values of Euclidean time xq. Calculation shows that even at the 
xq = 1/(2T) the tree level contribution is ~ 85% of the total value of the correlation function. Note also that large uj 
behaviour of the spectral function leads to a fast decrease of the correlation function C{xo) ^ l/a^o for small xq- For 
this reason the signal/noise ratio for the C(xo) is small at a:o ^ a and lattice measurement of the correlation function 
at Xq ~ 1/(2T) becomes computationally very expensive. 

In numerical simulation we use Wilson gauge action for the SU(2)-gluodynamics 


Sg=/3 




1 - -Tr{7^,^(x) 
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where U^^vix) is the product of the link variables along elementary rectangular (/r, u), which starts at x. 
For the tensor F^i, we use the clover discretization scheme: 

— Uu^^{x)). 


( 7 ) 
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FIG. 1: Renormalized correlation functions C{xo) as a function of Euclidean time xo measured at lattices 16 x 32® and 18 x 44® 


It causes no difficulties to build energy-momentum tensor having expression for the tensor F^^,. 

To calculate shear viscosity one should measure the correlation function Q. To carry out this measurement 
we use two-level algorithm described in M- This algorithm significantly improves the speed of the calculations. 
Note also that instead of the correlation function {Ti 2 {x)Ti 2 {y)) in this paper we measure the correlation function 
\[{Tii{x)Tii{y)) — {Tii{x)T 22 {y)))■ Both correlation functions are equal in the continuum limit[S]. 

It has become conventional to present the value of shear viscosity as a the ratio viscosity-to-entropy density 77 /s. 
For homogeneous systems the entropy density s can be expressed as s = , where e is the energy density and p is 

the pressure. These thermodynamic quantities were measured with the method described in m- 

Energy-momentum tensor in continuum theory is a set of Noether currents which are related to the transla¬ 
tion invariance of the action. In lattice formulation of field theory continuum translation invariance does not ex¬ 
ist and renormalization for energy-momentum tensor is required. For the correlation function considered in this 
paper the renormalization is multiplicative [18) . Renormalization factors depend on the discretization scheme. 
For instance, for the diagonal component of (when p = v) and the plaquette-based discretization of 

T/m = 77 ^ I — X) Trl7^^^(a;) -I- ^ TiUcy^y{x) j the renormalization factors are related to the anisotropy coef- 

ficients [1^1110]: ^ ^(piaq) — i_ — c,-), where and are defined in [T7] and for SU(2) 

are computed on the lattice in EH- 

Using the renormalization factors for the plaquette-based discretization of Tqo, we can find the renormaliza¬ 
tion factors for the clover discretization simply by fitting the vacuum expectation values of the renormalized Tqq: 

^{plaq) l^{p^^q)\^ _ ^(c/ot;) ^ 


Numerical results 

We carried out the measurements of the correlation function C{xo) at lattices 16x32® with (3 = 2.81 and 18x44® with 
P = 2.85. Both lattices correspond to deconfinement phase with the temperature T/Tc ~ 1.2 but different physical 
volumes and lattice spacings. In Figj^we plot renormalized correlation functions 0 as a function of Euclidean time 
xq for the lattices 16 x 32® and 18 x 44®. Two-level algorithm allowed us to reach accuracy ~ 3% for the lattice 
16 X 32® and ^ 10% for the lattice 18 x 44® at point Txq = 0.5. For the other points the accuracy is much better. 

Visually the set of values of the correlation functions measured at different lattices lies on one curve. The values 
for the correlation functions measured at close points in the region Txq ~ 0.5 are very close to each other. Note also 
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FIG. 2: The ratio piat{oj)/Pcont{^^) as a function of ua 


that the fit parameters (see below) of both correlation functions are very close to each other. These facts allow us to 
expect that finite lattice size and finite volume effects are not very important in our calculation. 

Unfortunately the accuracy of our results and number of points where the correlation function Q was measured 
do not allow us to use Maximal Entropy Method to determine the spectral function p{uj). For this reason we use 
physically motivated ansatz for the spectral function with unknown parameters which will be determined through the 
fitting procedure. Probably the simplest formula for the spectral function inspired by QCD sum rules [22] can be built 
if we join hydrodynamical behaviour at small frequencies with asymptotic freedom at large frequencies ^ 

Pi(w) = BT^ UJ d(uJo - w) + Apiat(uj) 0(uj - Wo). (8) 

In last formula piat(M>) is a tree level lattice expression for the spectral function calculated for the correlation function 
~ \{{Tii{x)Tii{y)) — {Tii{x)T 22 {y))) with clover discretization of the tensor at lattice with fixed Lt and Lg 
oo. Although the calculation of the piati>-^) can be easily performed using formulas from paper |23j . the resulting 
expression is very cumbersome. For this reason instead of the explicit expression for the piat, in Figj^we plot the 
ratio piati^)/Pconti^) for the lattice Lt = 16. 

The fit of the lattice data {xo/a > 2) with the formula (§ with the spectral function (|^ gives A = 0.723±0.002, B = 
0.079±0.016,wo/T = 7.5±0.5, xVdo/ ~ 1.4 for the lattice 16x323 and A = 0.703±0.003, B = 0.096±0.026, wq/T = 
8.3 ± 0.7, /dof ~ 0.7 for the lattice 18 x 44^. It is seen that ansatz ([^ for the spectral function fits the lattice 

data very well and the parameters of the fit for different lattices are in a reasonable agreement with each other. The 
resulting shear viscosity is p/s = 0.179±0.036 for the lattice 16 x 32^ and p/s = 0.217± 0.059 for the lattice 18 x 44^. 
In Figj^we plot the spectral function pi(w) as a function w for the lattice L* = 16. 

Now few comments are in order 

• In formula ^ we used piat{^^) instead of the continuum tree-level expression for the p(w) (see formula Q). We 
believe that this allows us to take into account discretization uncertainty and asymptotic freedom contribution 
more carefully. If one puts the continuum tree-level expression to formula ([^, the x^/do/ of the fits will be 
considerably enhanced: x^/do/ ~ 5.2 for the lattice 16 x 32^ and x?/dof ~ 4.5 for the lattice 18 x 44^. 

• As was noted above the asymptotic behaviour of the spectral function at large w is fixed p ^ w'^. However, the 
coefficient in front of the w^-behaviour is modified by higher order radiative corrections (see formula ([^). In 


Note that the frequency u) is measured in physical units. 
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FIG. 3: Spectral functions pi(uj) and P 2 (<^) as the function of tu/T for the lattice Lt = 16. 


order to take this effect into account in formula ® we introduced factor A. At this point our ansatz deviates 
from that used in papers [M E]. If we take A~^= 1, the ansatz Q will not be able to describe our data 
(x^/dof > 100 for both lattices). 

• The values of the fit parameters are physically well motivated. For instance, the value of the parameter A is 
smaller than unity, what agrees with the next-to-leading order result The value of the strong coupling 
constant at the threshold parameter wq (wq ~ 2.7 GeV in physical units jis q;s(wo) ~ 0.2 — 0.3. This allows us 
to expect that perturbative expression for the spectral function is applicable for to > ujq. 

Low frequency part of the spectral function Q is given by the first-order hydrodynamic expression ~ co. Comparison 
of the spectral functions of the energy-momentum tensor correlation functions obtained in N = 4 SYM [24] and the 
first-order hydrodynamic expressions allows us to expect that this approximation works well up to cu < ttT ~ 1 
GeV|25j. From the other side high frequency perturbative expression for the spectral function is fixed very accurately 
and it works well for w > wq = 2.7 GeV. The form of the spectral function in the region 1 GeV < w < 2.7 GeV is 
not clear. Formula Q continues the first-order hydrodynamic expression to the region 1 GeV < w < 2.7 GeV (see 
Fig. §. As the result there is rather large discontinuity in the spectral function (|^ at point ojq which for the lattice 
16 X 32^ is (pi(a;o -I- 0) — pi(a;o — 0))/pi{ujo — 0) ~ 6.7. The last fact allows us to state that the spectral function in 
form ([^ underestimates real spectral function in the region 1 GeV < w < 2.7 GeV and as a result the value of shear 
viscosity obtained in this fit is larger than its real value. 

As was noted a serious drawback of the spectral function p(w) is its discontinuity at point uj = ujq. In order to solve 
this problem we propose another form of the spectral function which preserves basic ideas of the spectral function Q 
and takes into account the property p{—oj) = —p(x) 

P 2 {uj) = BT^u) + Apiat{uj) tanh^ 

The fit of the lattice data {x^/a > 2) by the formula (§ with the spectral function © gives A = 0.723 ± 0.003, i? = 
0.039±0.012,a;o/T = 5.6±0.6, x^/dof ~ 1.2 for the lattice 16x323 and A = 0.705±0.004, B = 0.055 ± 0.022, wq/T = 
6.5 ± 0.9, x^/dof ~ 0.6 for the lattice 18 x 44^. It is seen that ansatz ([^ for the spectral function fits the lattice 
data very well and parameters from different lattices are in a reasonable agreement with each other. The resulting 
shear viscosity is 77 /s = 0.088 ± 0.027 for the lattice 16 x 32^ and 77 /s = 0.125 ± 0.050 for the lattice 18 x 44^. In 
Figj^ we plot the spectral function p 2 (w) as a function w for the lattice L* = 16. Notice that the function Q is 
larger than ([^ in the region 1 GeV < a; < 2.7 GeV, as the result the value of the of shear viscosity becomes smaller. 







6 


Note also that the value of shear viscosity obtained at the lattice 16 x 32^ is very close to the N = 4 SYM prediction 
ry/s = l/ 47 r ~ 0.080. 

In order to investigate the systematic errors, caused by the usage of the ansatz ( 8 ),(9) we studied other fit functions. 
Instead of the tanh^(a;/a;o) in formula (|^ one can use, for example, any power of this function ^ tanh^^(a;/a;o) or 
linear combination ^ A*, tanh^^(a;/a;o). ^ We found, that in these cases the fits are also good and the resulting 

viscosities are greater than that for the spectral function and smaller than that for the spectral function ([^. 

As the result of this paper we take the value 

- =0.134 ± 0.034 ± 0.046. (10) 

s 

The central value is the average between the values of shear viscosity obtained with the spectral functions in the form 
([^ and from the fit of lattice data measured at the lattice 16 x 32^. The first uncertainty is due to the fitting 
procedure which is typically 25%. The second uncertainty is due to the unknown model of the spectral function. We 
estimated this uncertainty assuming that the central value of shear viscosity can vary between central values of the 
fits ([^ and (|^. Actually there are a lot of different sources of uncertainties of our result but they are much smaller 
than statistical and spectral function model uncertainties. 


Conclusion 

In this paper we studied the energy-momentum tensor correlation function (Ti 2 ( 0 )Ti 2 (cc)) for 5't7(2)-gluodynamics 
using lattice simulation of QCD. We carried out the measurements of this correlation function at lattices 16 x 32^ 
with P = 2.81 and 18 x 44^ with /3 = 2.85. Both lattices correspond to deconfinement phase with the temperature 
T/Tc ~ 1.2 but different physical volumes and lattice spacings. In order to enhance the accuracy of the calculation we 
used two-level algorithm which allowed us to reach accuracy ^ 3% for the lattice 16 x 32^ and ^ 10% for the lattice 
18 X 44^ at point xq = Xj^T. For the other points the accuracy is much better. 

Using lattice data for the correlation function we tried to study the spectral function. It was shown that physically 
motivated anzatz which joins the first-order hydrodynamical behaviour at small frequencies with asymptotic freedom 
at large frequencies fits our data very well for both lattices. We also studied other forms of the spectral functions 
which interpolate between hydrodynamics and asymptotic freedom. 

The results of the study of spectral functions allowed us to estimate the ratio of shear viscosity to the entropy density 
77 /s = 0.134 ± 0.057. This value is in agreement with our previous finding 77 /s = 0.111 ± 0.032 [12] obtained at small 
lattice with the estimation of statistical uncertainty only. Note also that within the uncertainty of the calculation the 
value obtained in this paper agrees with the value measured in S't7(3)-gluodynamics: 77 /s = 0.102 ± 0.056 measured 
at temperature T/Tc = 1.24|10j. 

The values of the ratio 77 /s for the SU{2) and S'?7(3)-gluodynamics are very close to the prediction of the N = 
4 SYM theory at strong coupling: 77 /s = l/ 47 r. So one can state that both theories belong to strongly correlated 
systems and probably the closeness of the ratio 77 /s to the ratio in N = 4 SYM theory is the property on non-abelian 
gauge theories. 
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